NPS ARCHIVE 
1963 

MAULDIN, 3 . 

(graduate Aeronautical Laboratories 




Firestone Flight Sciences Laboratory 



Guggenheim Aeronautical Laboratory 



'an Laboratory of Fluid Mechanics and Jet Propulsion 

Thesis 

M384 



Pasadena 



I ART I. ANALYSIS G P A MULTI-STAGE AXIAL CGMi- RESSCR 
WITH HIGH REACTION BLADING 



PAR.T U. A DESIGN STUDY GF A MULTI-STAGE AXIAL 
COMPRESSOR -WITH BLADING OF HIGH ASPECT RATIO 



Thesis by- 
James H. Mauldin 

// 

Lieutenant, United States Navy 



In Partial Fulfillment of the Requirements 



For the Decree of 



Aeronautical Engineer 



California Institute of Technology 



Pasadena, California 



^rSeschoou 
TEREV CA 93943-5101 



1963 






I v_i- V ^ 

J . /I • 



7T 



ACKNGV- JUDGMENTS 



The author wishes to ejcpress his sincere appreciation to D 
V. D. Rannie, who originally suggested the investigations in this 
thesis. Eis many hours of encouragement, guidance, and helpful 
criticism were invaluable during the course of this work. 

Thanks are also due to Mrs. Roberta Duffy for typing the 



manuscript. 



PART I. 



ANALYSIS OF A MULTI-STAGE AXIAL COMPRESSOR 
WITH HIGH REACTION BLADING 



ABSTRACT FOR i-ART I 



An analysis of a multi-3tage axial compressor with high reac- 
tion blading is carried out. The usual methods of linearizing the equa- 
tions for flow through the compressor fail for this type of blading. A 
numerical solution of the non-linear equations is worked through for 
off-design operation. Flow velocities through the first six stages are 
calculated for the off-design flow rates in order to insure that the flow 
does indeed come to a steady* state repeating pattern. Entering flow 
angles for both rotor and stator are calculated for this repeating flow 
condition. The off-design incidence angles for the stators at low flow 
rates indicate a possible deterioration of efficiency. 



SUMMARY FOR PART I 



The pressure rise per stage in an axial flow compressor is pro- 
portional to the square of the incidence velocity relative to the blades. 
Almost all axial compressors for gas turbines were designed with the 
whirl in the same direction a3 the rotor blade speed. This choice of 
design flow gave high efficiencies when the blades were designed prop- 
erly. Until the introduction of transonic blading, a high efficiency was 
obtained only when the Mach number relative to the blade was less than 
unity -- usually in the neighborhood of 0. 75. Relative velocities of this 
magnitude did not require whirl opposite to the direction of rotor blade 
speed and, in fact, the matching with a turbine was easier with whirl in 
the same direction as the rotor speed. 

An axial compressor for a gas of low molecular weight, for in- 
stance helium, can be designed with blading similar to the axial compres- 
sor for air. However, the pressure rise per stage is proportional to 
the square of the Mach number of the relative velocity. The rotor blade 
rotational speed cannot be increased because of strength limitations, and 
hence the square of the relative Mach number and therefore the pressure 
ratio are decreased by a factor approximately proportional to the ratio 
of molecular weights. Thus, the number of stages required for a com- 
pressor using helium would be approximately seven times that for air 
for the same overall pressure. Clearly, higher relative velocities are 
required to decrease the number of stages for a gas of lower molecular 
weight. An obvious method of doing this is to introduce the whirl in the 
opposite direction to rotor blade speed. 

A vortex flow pattern was chosen as the design condition: then 
the flow is truly two-dimensional; that is, there is no radial shift of the 
streamlines as the fluid proceeds downstream through the compressor. 



This considerably simplifies the analysis of the flow at the design flow 
rate. Off-design calculations must still be carried out in order to in- 
vestigate the compressor performance and to indicate how rapidly losses 
might increase as the compressor flow rate increases or decreases from 
the design value. For conventional multi-stage compressors with rela- 
tively low ideal pressure coefficients, this analysis is usually carried 
out by linearizing the governing equations for three-dimensional flow 
through the compressor, and is quite satisfactory. However, with whirl 
against the rotor blades and its correspondingly higher pressure coef- 
ficient, this type of linearization fails completely. It is necessary to al- 
low for the streamline shift to avoid singularities in the solution. Thus, 
the full set of non-linear equations, talcing into account the radial shift 
of the streamlines, must be used. In order to calculate the repeating 
flow pattern characteristic of multi-stage compressors, it was necessary 
to numerically determine the flow pattern through the first four or six 
stages. 

Although detailed calculations of efficiency are scarcely possible, 
the results of the analysis can be used to give indications of probable 
losses. The degree of turning required by the blades appears to be with- 
in permissible limits. The off-design incidence angles for the rotors 
deviate very little from the design values. The stator incidence angles 
are more severe, and might produce a loss in efficiency over part of the 



stator blades. 



I. INTRODUCTION 



Conventional subsonic axial-flow compressors designed to U3e 
air as a working fluid are not necessarily aiitable for gases of lower 
molecular weight. For a particular closed cycle system, where a 
working fluid other than air might be desirable, the number of stages 
required will vary greatly with the choice of gas. The number of stages 
is directly proportional to the specific heat at constant pressure and is 
inversely proportional to the blade speed and the change in tangential 
velocity of the gas through a stage. The blade speed is fixed by the 
strength of the blade material and is generally held as high as possible. 
With present-day air compressors, there i3 a maximum amount of 
tangential velocity change that can be accomplished efficiently. Thus, 
for a typical compressor designed to use air, the number of stages re- 
quired would be increased approximately five times if a gas such as 
helium were used and the same temperature ratio across the compres- 
sor were to be maintained. This report describes a blade design that 
will hopefully reduce this high number of stages required for helium to 
approximately one-half the number presently used. 
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II. THEORY 

Unlike conventional multi-stage compressors, the design of the 
stator blade is such as to give an accelerating flow through the stator 
rows. The rotor bladee then turn the flow back to the axial direction. 
This type of design has been used in some single-stage fans. The nota- 
tion to be used in this report is illustrated in Figure 1. 




Figure 1 . 

The following list of symbols will be used throughout this report: 
c a = axial component of velocity 
c u = tangential component of velocity 
c = absolute velocity 
w = velocity relative to the rotor 
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vv u = tangential velocity relative to the rotor 
•y = angle between axial direction and absolute velocities 

p = angle between axial direction and relative velocities 

R = tip radius 

u = tip speed 

cij = angular velocity of rotor 

r = radius measured from compressor centerline 
p = static pressure 

p = density (assumed constant throughout) 

Subscripts refer to stations as shown in Figure 1. 

In "two-dimensional" flow through an axial compressor, the 
streamlines are assumed to lie on cylindrical surfaces with constant 
radii. The pressure gradients in the radial direction are ignored, al- 
though it is shown later that at design for vortex blading, the flow is 
two-dimensional and the pressure gradients are properly accounted for. 
Since an important difficulty in applying two-dimensional flow theory to 
this design arises, we will discuss it briefly here. 

The blades are of the free vortex design. Thus 

c « 1/r . 

U 1 

A dimensionless radius is defined by 

£ = r/R . 

Therefore, the tangential component o£ velocity through the stator is 
taken as 

c = ua/£ , 

U 1 

where a is an arbitrary constant to be determined. 

The pressure drop across a stator row is given by 



( 1 ) 
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(2) 



Since r is a constant along a streamline, the Bexmoulli equation for 
rotating axes gives, for the pressure rise across a rotor row, 



p :- " p i , , 2 

— = " < w i - w 2 ) • 

The pre3suierise across a stage is then 



( 3 ) 



Po - P 



o i/2 2 2 2 1 

- = rlwj -w 2 +c a -cj ). 



( 4 ) 



At design, c^ = 0 , and from the velocity diagrams we obtain 



W = U -Jr- + u£ 



u 



f. 



W = U.: 

u 2 

2 2 2 2 2 a „ .2 

w, = c +w = c + u (r- + g) 

a 5 



1 



u 



1 



2 .. 2 



L £. C L , Cg. 

v/_ = c +w = c +u£ 

/ rs a ^ 



U- 



2 . a .2 



= c„ + c = C + u ( F ) 



u 



After substitution of these relations into ( 4 ), we obtain, for the pres- 
sure rise across the stage, 



'2 *o 2 

= u a . 



A work coefficient is usually defined as 

A? 



* = • 

yp uu R 



( 5 ) 



(6) 



Thuo, when the compressor is operating at the design point, 

* = 2a . 

However, a is not unrestricted. It has been found experimentally 



desirable (Ref. 1 ) to have 



-5- 



W 1 w 2 



w , 



< P . 



(7) 



where p is between 0. 5 and 0. 6 . Defining a non- dimensional flow 
coefficient as 

9 = c 3 / u • (8) 



this restriction can be written 

2 ^ 2 
(» + 2|J> A 

4p/i-p P/p 



(9) 



This condition must be satisfied over the entire length of the blade. 

Since the restriction is most severe at the hub, values of t and <$> 
evaluated there will be conservative over the rest of the blade. An ex- 
isting experimental set-up prescribed the hub ratio = 0. 6 . There is 
another restriction on the velocities that is not easily described in analy- 
tical form. This is the amount of turning that can be accomplished ef- 
ficiently in the blade rows. Due to lower blade speed at the hub and the 
adverse pressure gradient in the rotor, this condition also is most 
severe at the rotor hub. p was chosen as 0. 6, and values of the turning 
angle were obtained for different values of ’ . The quantities ijr = 0. 70 
and <j> = 0. 58 were subsequently chosen. The value of i' can be raised 
substantially for a compressor with hub ratio in the range 0. 7 to 0. 75. 

So far, no difficulty with the two- dimensional theory has arisen. 
This is not unexpected since, as mentioned before, vortex blading is 
exactly two dimensional at the design flow rate. However, we have yet 
to consider the off-design conditions where a radial component of ve- 
locity does exist. At design, the work coefficient was constant over the 
blade height. This is desirable in order to cut down losses due to mix- 



- 6 - 



ing at the end of the compressor, and actually is necessary on all stages 
after the first few because a piling up of pressure at some radius is 
impossible for a very large number of stages. Thus, this restriction 
will be applied at off-design flow rates also. Then 



Ap 1 , 2 2, 2 2. . . 

■- -TT = -TTT (W 1 * w 2 +c 2 ‘ C 1 ) = constant. 

yp «! R Ui R 



( 10 ) 



Remembering that c^ 4 0 for off-design conditions, then 



w = uS, + c 



u 



1 






w = u£ + c , 
u, ^ u> 

andlhe work coefficient becomes 

if' = 2 — (c -c ) 
u U 1 u 2 



or 



t’j = 2£f£ + <j>(tanyj - tanf^)] • (11) 

A first approximation can be made that tany ^ and tan j?> 2 , correspond- 
ing to the leaving angles of flow, are independent of the inlet angles and 
hence are independent of flow rate. If this is done, the local value of 
at the root is 



= 0 . 6 ) 



0.72 - -j> 

1. 2(tan s Z, - tany j) 



( 12 ) 



Since under our approximation ( tanO^ - tany^) is constant and is great- 
er than zero, we see that for * > 0. 72 the flow coefficient at the root 
becomes negative and is meaningless. If tan'l^ and tany ^ are allowed 
to change with flow speed and this change calculated by cascade theory, 
the singularity in is not relieved appreciably. Since values of the 
ideal pressure coefficient greater than 0. 72 would be required in starting. 
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for instance, the two-c&nensional analysis of the flow is unsatisfactory. 
It must be mentioned that if * could be pushed as high as 0. 72 in a 
standard compressor, there would be the same difficulty. Normally, 
axial-flow air compressors have an ideal pressure coefficient of be- 
tween 0. 4 and 0. 5 . It will be shown that a three-dimensional flow 
analysis of the compressor does not contain this singularity in § . 

This difficulty is mathematical rather than physical, and appears to be 
the result of an invalid linearization of the three-dimensional equations 
of motion. There seems to be no way of calculating the ideal flow pat- 
terns except by a numerical procedure with a radial shift of the stream- 
lines taken into account. 

The assumptions made in the three-dimensional flow analysis 



are: 



(1) incompressible perfect fluid; 

(2) axisymmetric flow; 

(3) equations are written for conditions far upstream and far 
downstream of a blade row. 

The notation used is the same as in the two-dimensional analysis 
except that the odd subscripts refer to stations far downstream of stator 
rows and the even subscripts refer to stations far downstream of rotor 
rows (see Figure 2). 

The Bernoulli equation for flow through a stator row is 



2 P 



2n+l . 2 

— — + c 
P u 



+ c 



2n+l 



2n+l 



2? 2n 2 2 

T C + C 

P u-,„ a 



(13) 



2n 



2n 



The Bernoulli equation for flow through a rotor row is 
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Station 2n+l 






Station 2n 




dr Zn+l 
















dr 7 

2n 








blade 

row 









r 2n+l 



Figure 2. 



2p 



2n . 2 2 2 2 

+ w + c - id r 0 

P u, a-,„ 2n 



2p 



2n 



2n 



2n+l . 2 

+ w 

p u 



2 22 .... 
+ c a ' r 2n+H 14) 

2n+l a 2n+l ^ n+1 



The equations for the pressure gradient in the radial direction are 

. 2 2 

dp 2n d P?--' Cu - 



ht: 



= p 



2n 



and 



2n+l 



Hr! 



= P 



‘2n+l 



(15) 



'2n 2n 2n+l 2n+l 

By using (15) to eliminate the pressure terms, equations (13) and (14) 



become, respectively, 
2 



u 



2n 



2n + 1 



73 ^ (c 4 +c 4 ) j dr2 " 



u 



2n+l , 1 



2n+l 



73F 



2n+l 



( c u 



+C 

2n+l a 2n+l'J 



dr 



2n+l ’ 



(16) 



u 



2n + 1 



t zxi 



d 

m ■■ 

*-> 
2n 



2 2 „ 2 2 v 

+c -to r ) 

cn a 2n 2n ' 



dr 



2n 



r.C 



u. 



2n+l A 1 



2n+l 



(w ; 

\ u 



+c 

2n+l a 2n+l 



2 2 
-tu r 



2n+l/ 



dr 



2n+l 



(17) 
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Considering a cylindrical stream tube, the continuity equation can be 

written in the differential form 

c r, dr., = c r, ,.dr, , , . 

a^ n 2n 2n a 2n+l 2n+1 2n+l 

Within the scope of the assumptions previously mentioned, these equa- 
tions are quite general. It is now assumed that the leaving angles of 
both stator and rotor are constant with respect to flow rate and give 
vortex flow at the design flow rate. Thus, evaluating the leaving angles 
at the design condition, indicated by subscript d , we obtain 



u 



* C 



2n+l 



2n+l 



tanvj = - 



R 

2n+l 



(18) 



or, in the non-dimensional notation, 

♦ 



tany . = 



a 



W. 



d 1 

in 



•2n+l v d * 2n+l 
Similarly, we obtain for the rotor leaving angle 

*2n 



tanjB^ = 



(19) 



( 20 ) 



With the leaving angles so specified and tlie following definitions, 

and b » , 



a = 



d 



the equation for the flow through the stator becomes 



( 1+ ?H 



2 ‘ £: 



2n+l _ ^2n+l 



^2n+l 7 ^ 2n+1 



’2n 



L d *2n 



( 



i 2^2 

b +£ 



_i_ 



^2n) ' 



2| 



2n 



7~7T 



b +£ 



2n 



(Jl. 

V«2„ 



J 



The equation for flow through a rotor row becomes 



- 10 - 



r b 2 +| 2 






2n 



2n 

2 ^2n 



2| 



2n 






e 



2n 



*5,-1 (L“ a 2 +| 2 n -l < ^ 2n ~ 1 

The continuity relations are 

^2n- 1^ 2n- l d ^ 2n- 1 = < * > 2n^2n d ^2n = ^2n+1^2n+l d ^2n+l ’ 



1 + 



a.£ 



B — 
2 

2n-l 



1 + 



^n-d 



d *2n-l 



( 22 ) 



(23) 



The two pairs of non-linear equations (21) and (22), and (23) can be 
solved approximately for <j> m as a function of for m = 0, 1, 2, . . . 

in succession. 

Notice the particularly simple solution 4> m = ^ for all values 
of m . It follows that there is no streamline shift and no radial velocity 
in this case, and the analysis of the flow reduces to the two-dimensional 
one previously considered. It is largely because of this analytical sim- 
plification that the free vortex design has been used so frequently. How- 
ever, for off-design flow rates, no svich simple solution exists, even 
for the free vortex design. 

To illustrate the method of solution to these equations, let us 
consider an off-design constant axial flow approaching the first sta- 
tor. In this case, care must be taken in applying equation (21), since 
the assumption of a whirl velocity is incorporated in the right-hand side 
of this relation. For this incoming flow, equation (13) becomes 



2p 



1 . 2 2 
- + c + c 
a, u. 



2p 



— + c = constant, 
a 

o 



Since dp^/dr^ = 0 , equation (16) becomes 



u. 



Hr 



(c + c ) = 0 . 
U 1 a l 
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Substituting from equation (19), the velocity diagram relations, and the 

definition of a , the equation corresponding to equation (21) is 

2. d*. 



( 



1 + = o 



if' 3 *! 



or 



6 j - constant. 

This constant is then determined by integrating the continuity relation 

from the hub to the tip. 

1 1 

J*l*l d *l = /W*0 
1(1 - 6*1 = « 0 1(1 - 6 *) 



*i = 



Equations (21) and (22) are applicable to all of the other blade 
rows now. The relation for <{>2 becomes, using equation (22), 



d 



,2 . 2 

b +£. 



2|- 



= 0 . 



Integrating, we obtain 

*2 * b 





where K is determined as before by integrating the continuity expres- 
sion. 



The flow coefficients for the following blade rows are obtained in 
the same maimer; however, the integration for these must be done nu- 
merically. Also, another mathematical difficulty arises in the solution 
for (j> 3 and the remaining stages. The existence of both independent 
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variables on the right-hand side of both { 21 } and ( 22 ) will require an 
iteration process. The procedure is to first let ^ * £gn solve 

for $2^1+1 • The continuity relation can now be integrated from to 
£ , and a relation between ^n+l ^Zn °k ta toed. Substitution of 

this relation into the differential equation for §2x1+1 y^ e ^ s a refined 
value for ^ 2 n+l * the correction to the flow coefficient is 

not large, and the iteration process is generally not necessary. 
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in. APPLICATION OF THE THREE-DIMENSIONAL THEORY 



The calculation of the flow field as it progresses down the com- 
pressor was carried out for two off-design conditions. The blades were 
chosen such that <j>^ = 0. 57751 and = 0. 70000 . The average flow 
rate for the off-design conditions was taken as <j> = . 46000 and $ = 

. 69000 . Figures 3 through 6 show the variation of the flow coefficient 
versus dimensionless radius at various stations in the compressor. 

The flow field at lower than design flow rate is converging fairly rapidly 
to a steady-state condition, while at the higher flow rate the convergence 
does not appear to be as good. 

The error in choosing ^n+l = ^2n as a ^ rst approximation is 
shown in Figure 7. The calculation was carried out for the fourth rotor 
row. An error of only approximately 2 percent of <|>g at the mid-height 
of the blade is introduced by the first approximation. It is felt that real 
fluid effects will be greater than this amount. 

The local values of * for the two off-design flow rates are shown 
in Figure 8. Notice that is very nearly independent of the radius, 
agreeing with our assumptions leading to equation (10). 

In standard compressors there exists a perturbation solution to 
the three-dimensional equations. There it is assumed that the flow 
from a stator is approximately equal to the flow ^ from a rotor. Thus 
the procedure is to find a solution for which ^ « (4> 2 + ! i > j) • To 

the zeroth order, this results in a linear compressor performance 
curve in the form 



* = A - B4> . 



(24) 
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The two-dimensional theory also produced an equation of this type, 
since in that case was taken identically equal to 4*2 * 

To examine the possibilities of this approach for the proposed 
helium compressor, consider the approximate performance equation 



t? = 2[£~ + £cj>jtanyj - j^tan^] . (25) 

For free vortex blading, equation (25) becomes 

= 2 757 + • < 26 > 

If we now assume 

+ ^ + < 27a > 

and 

$2 = 7<V<>2> - i(4> 1 -4‘ 2 ) n <> 0 - A 4> . <27b) 

then equation (26) can be written 

* * 2 s 2 + ^[V 2 * 2 } + x: l >d +2 S 2 i • < 28 > 

v d v d 



Typical values of and A<|> for a conventional compressor allow the 

term containing A<f> to be ignored. However, for the compressor pro- 

2 

posed in this report, t ^ is approaching 2£ at the hub, and the term 
containing Acj> can no longer be neglected. Figure 9 shows the per- 
formance curves of both types of compressors with and without the in- 
clusion of the term. The inadequacy of the two-dimensional 

solution in the case of the helium compressor is obvious. To the scale 
of the graph, no difference can be observed for the conventional com- 
pressor. 

Figure 10 shows the total pressure variation along the blades 
after the flow has essentially settled into a steady pattern. For lower 
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than design flow rates, more work is being done by the tips, while at 
higher flow rates, the roots give the greater rise. The cumulative dif- 
ference after six stages is still less than half the design pressure rise 
for a single stage. 

It was mentioned previously that it is desirable to have 

2 2 
w , - w _ 

— < 0 . 6 . 

v 'l 

This quantity is plotted in Figure 11. Possible difficulties are indicated 
at lower than design flow rates. A higher hub ratio would eliminate this 
difficulty. The corresponding ratio for the stator is plotted in Figure 
12. Since the stators are in a field of decreasing pressure, no diffi- 
culties are expected here. 

The off-design entering angles for the stator and rotor are plot- 
ted in Figures 13 and 14, respectively. The stators are very critical 
and will require care in their design. Sharp leading edges of the sta- 
tors would probably cause high losses. The off-design flow incidence 
angle deviations for the rotors appear satisfactory. 
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IV. CONCLUDING REMARKS 

It has been demonstrated that high reaction blading can be used 
in multi-stage axial compressors to produce a significant increase in 
pressure rise per stage over conventional blading. This type of blading 
could be used in a closed-cycle helium compressor and greatly reduce 
the total number of stages required. The analysis also indicates that 
a hub ratio greater than 0. 6 would possibly be a more efficient design, 
since stator incidence flow angle deviations become quite large near the 
hub in the example shown. 

It was also shown that the usual linearization of the three- 
dimensional equations of 'flow is not valid for this type of blading. How- 
ever, a numerical solution allowing for the shift of the streamlines does 
converge to a steady, repeating flow pattern after the first five or six 
stages for flows as much as twenty per cent off design. 
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Figure 3. Stator Axial Velocity Coefficients, $ = . 46000 . 
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Figure 4. Rotor Axial Velocity Coefficients, $ = . 46000 . 
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Figure 5. Stator Axial Velocity Coefficients, <£> = . 69000 . 
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Figure 6. Rotor Axial Velocity Coefficients, <J> = . 69000 . 
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Figure 7. Streamline Radii Upstream and Downstream of a Rotor. 
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Figure 8. Local Value of $ Versus Non-Dimensional Radius. 
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Figure 9. Effect of Two-Dimensional Linearization on Performance Curves. 
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Figure 10. Total Pressure Variation with Radius. 
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Figure 11. Kinetic Energy Conversion Through the Rotor Versus 

Radius. 
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Figure 12. Kinetic Energy Conversion Through the Stator Versus 

Radius. 
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Figure 13. Stator Flow Entering Angles. 
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Figure 14. Rotor Flow Engering Angles. 
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Figure 14 (continued). Rotor Flow Entering Angles. 
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PART II. A DESIGN STUDY OF A MULTI-STAGE AXIAL 
COMPRESSOR WITH BLADING OF HIGH ASPECT RATIO 



- 32 - 



ABSTRACT FOR PART II 

A method for designing blades of high aspect ratio in a multi- 
stage axial compressor is investigated. This method requires a cal- 
culation of the flow induced by all blade rows. The method is illustrated 
for the limiting condition of hub ratio equal to unity and then is carried 
out for the more realistic case of finite hub ratio. 

An example of a blade design for a particular flow is carried 
out. The results are compared with those of a previous design for the 
same flow conditions but based on a theory applicable for blading of low 
aspect ratio. 
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SUMMARY FOR PART II 

The flow through most present-day axial flow compressors is 
frequently analyzed by applying a three-dimensional theory that is based 
on the assumption of blading with low aspect ratio. Even though existing 
compressors do not have what could be called low aspect ratio blades, 
the theory predictions and experimental results show reasonable agree- 
ment. 

Where weight considerations are of primary concern, the possi- 
bility of reducing the size of a compressor by shortening its length is 
attractive. There is no known reason why the blade rows cannot be held 
to their same height but shortened in chord. There is, however, a dif- 
ficulty in applying the aforementioned analysis. There, it is assumed 
that equilibrium flow patterns exist at the leading and trailing edges of 
the blades. Tlis is identical to assuming low aspect ratio for the blades, 
since they must be wide enough to allow all changes in the flow to occur 
entirely within the blade row itself. Several extensions of the theory 
have been developed to take account of the non-equilibrium flow near the 
blade rows. Some of these analyses assumed small vorticity throughout 
the flow field and also a small change in vorticity through the blade rows. 
Another theory relaxed the restriction of small vorticity throughout the 
flow field. 

The theory presented in this thesis is an extension of the latter 
approach, and is based on the assumptions of closely spaced blades of 
high aspect ratio. The equations of motion are linearized by a perturba- 
tion scheme that assumes small changes in perturbation velocities across 
the blade row3. To illustrate the method, the calculation is first 



- 34 - 



reetricted to the limiting case of hub ratio close to unity. Then the 
more involved problem for any hub ratio. is investigated. An example 
is worked for blades of aspect ratio about three. Here, there appeared 
little difference between the theory applicable to blades of low aspect 
ratio and the theory of blades of low aspect ratio. 
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I. INTRODUCTION 

The theory of three-dimensional flow in the axial compressor 
described in Part I of this thesis is concerned only with flow conditions 
far upstream and far downstream of a blade row. This theory would be 
applicable to blades of very low aspect ratio. The difference between 
these velocity profiles and the actual profiles near the leading and trail- 
ing edges of the blades was assumed to be negligible. It was on this 
basis that the entering flow angles for the blades were determined. If 
the aspect ratio is not very low, it is possible to obtain more exact in- 
formation about the flow pattern near the blade rows by making use of 
some of Marble's results (Ref. 1). This has been done by Bowen, 
Sabersky, and Rannie in Ref. 2. The method carried out by these in- 
vestigators forms the basis of the theory to be discussed in this report. 
The approach is to consider the disturbances due to rows of blades of 
high aspect ratio placed very close together. Here, the flow in any par- 
ticular blade row is influenced by neighboring blade rows. In the three- 
dimensional theory of Part I, the changes in the flow patte rn were as- 
sumed to occur entirely within the blade row itself and one blade row did 
not affect its neighbor. 

The analysis begins with the calculation of the disturbances due 
to an infinitesimally thin disk. These disturbances can then be super- 
imposed to form a blade with finite axial extent. The effect of many 
blade row3 can be determined by the superposition of the disturbances 
due to the finite row. In this way, the flow pattern for a multi-stage 

axial turbomachine can be constructed. For a repeating flow pattern, 
the summation of the effects of all blade rows is possible. 
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II. TWO-DIMENSIONAL THEORY 

A. Equations and Boundary Conditions 

Before proceeding directly into the three-dimensional theory of 
high aspect ratio blading, it is informative to investigate first the limit- 
ing case of high hub ratio. This problem is simpler because the radial 
pressure gradients vanish far from the blade. The analysis is based 
upon the disturbances due to a disk at x = 0 (see Figure 1), and the 
equations of motion are written subject to several assumed conditions. 
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Figure I 

These assumptions are: (1) there are no circumferential variations; 

9 

that is, ( ) = 0 . This is equivalent to having a disk with an infinite 
number of blades. (2) The flow is incompressible, the density p being 
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hcld constant. The continuity and momentum equations then become 



9u 3v _ fl 

3S + 3y " 0 




(1) 


Ou . 


Su 


1 Op 


(2) 


+v 3y = 


* 7 ix 


3v 


Ov 


1 Op 


(3) 


c 

ft- 

-f- 

< 

ft 

II 


‘pSy 


Ow 


8w 




(4) 


U -5— + V 75 = 

Ox Oy 


0 



where u , v , and w are the gas velocity components in the x , y , and 
z directions, and p is the static pressure. 

The boundary conditions to be imposed on this system are as 

follows: 

(1) v = 0 at y=0 and y = b for all x , since the walls of the 
compressor are streamlines. 

(2) u is continuous at x = 0 from continuity through the actu- 
ating disk. 

(3) v is continuous at x = 0 due to the component of force 
being zero in the y direction. 

(4) w is allowed to change discontinuously across x = 0 . 

(5) At x = - oo , 

u = u(y) - eu o (y), 
v = 0 , 

w = w(y) - fw Q (y) , 
p = p - cp Q = constant. 

(6) At x = -f oo , 

u = u(y) + ■' u Q (y ) . 
v = 0 , 
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w = w(y) + r\v Q (y) , 
p x p + cp o = constant. 

These last two boundary conditions are written in a form so as to con- 
form with the following perturbation scheme. Far upstream and far 
downstream of the disk, the flow is independent of x j any disturbance 
due to the disk having been smoothed out. Close to the disk, the flow is 
disturbed by a small amount, and the following perturbations are as- 
sumed, where the symbol c is carried only to denote the assumed 
small quantities. 

u x u(y)T cu Q (y) + cu^x.y) 

v x s v l (x t y) 

(5) 

w = w(y) + ew Q (y) + cw^x, y) 

p X p? rp Q + epj(x, y) 

Notice that p and p Q are assumed independent of either x or y and 
that the mean value of v is zero. The top signs are for - oo < x < 0 
and the bottom signs are for 0 < x < + co . 

Introducing equation (5) into the equation of motion yields the 
linear set 
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B. Solution of the Equations for the Actuating Disk 

The first three of the preceding equations (6, 7, and 8) can be 



combined to give the differential equation for v , as 

V 1 

"T 



a 2 v 



. 9 2 v. . n 2_ 

1 , 1 1 0 u n 

^ + T " ~ — 2 v i = 0 
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( 10 ) 



with the boundary conditions that v = 0 at y = 0 and y = b . A solu- 



tion to equation (10) that satisfies these boundary conditions is 
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+ X rr 
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where 



n=l 



U = 



u dy 



0 



and Y R (y) satisfies the differential equation 



d 2 Y 



n 



dy* 



u J 



Y = 0 
n 



with boundary conditions 



Y n<°> = V b > = 0 ‘ 



Equation (6) then yields 



co 



u, = T U Y A ~ e T Xn * Y* (y) 
A n \ n 

i n 

n=l 



( 12 ) 



where the prime denotes differentiation with respect to y . Notice that 
the boundary conditions on u are satisfied at x = - oo and x = + oo . 
rrom the continuity of u across the disk, the following relation is 
obtained 

oo 



u D (y) = - u £ A n ± Y^y) . 



n= 1 



(13) 
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Slnce u and u q are to be chosen and the shape of the blades is to be 
determined, it is possible at this point to choose = 0 for n > 1 . 
Also, if u is chosen as a linear function of y , then 




and 



Note that it is not necessary that these particular choices be made. It 
is done here only to reduce the amount of algebra in the rest of the analy- 
sis. It will be pointed out later that these choices can indeed represent 
the flow through a conventional compressor quite well. Therefore, the 
perturbation velocities become 
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Using equation (9) » Wj is determined as 
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and with the use of equation (7), 
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There are several requirements that we can now impose. First, 
we can require the total pressure rise across the disk to be a constant. 
This will insure that the flow pattern will repeat itself. This is a similar 
condition to that set upon the work coefficient in Part I of this paper. To 
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calculatc the change in total pressure, one needs the change in whirl 
and the change in static pressure. The change in wnirl is 
Aw s w(+o, y) - vv(-o, y) , 
and with the U3e of equations (5) and (14) becomes 



Aw = 2r\v -2c — 

o — 
u 



dw 

w 






u (r) dr 
o' ' 



The change in static pressure is 

s ^ <+©» y) - ^ (-o» y) . 

and using equations (5) and (15), becomes 
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Since u and v are continuous through the disk, the change in total 
pressure is 

Ap r 



+ w Aw 

P P 



or 
Ap, 



L » 2e ^2 + 2c (u u o + w w o ) + 2e (|j~ + ^ j 



du , w dw 
u 



( 16 ) 



and v/e require this to be constant. Thus, we now have the four unknown 
functions of y : u , u , w , and w , and two constants p and p Q 
with one equation (16) relating them. Also, since we are interested only 
in the pressure rise in a compressor, v/e are not immensely concerned 
with the mean pressure , p . The pressure change across the disk is 
?.p Q , and this is the more important quantity. 

Another condition that v/e can require is that the total pressure 
far upstream and far downstream be constant from root to tip. Note 
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that this ia not a necessary condition, but is a common choice in com- 
pressor design and gives two more relations connecting our four un- 
known functions of y . For the total pressure to be constant we must 
require that 

p 2 Z 

— + [(u-cu ) + (w-ew ) ] = constant 

p * o o 

and (17) 



+ £[(u+-u o ) 2 + (w+rv/ o ) 2 ] 



constant 



and since p/p is constant far upstream and far downstream of the 
dial?, we derive the following two equations by adding and subtracting 
equations (17) and keeping only first order in e : 



-2 t -2 
u + w = Kj 



(! 8 ) 

U U + W W a K, 

O O C 

where Kj and K-, are undetermined constants. 

The second condition of equations (18) is now actually the same 
as equation (16). If the disk were to act as a stator, then (16) would 
reduce to 
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(19) 



Since p /p is a constant, and from (18) uu + w w = K- , the last 
term of equation (19) must be a constant. Since the integral is obviously 
zero at y = 0 , then this constant must be zero. Thus, with the assump- 
tion of constant total pressure, we have four unknown functions of y 
and two conditions upon them. Hence, for example, we could pick u , 
u q , and p Q , and the flow will be completely determined. Once this 
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has been dene, the shape of the blades can be found. 

To derive the effect on the flow due to a blade of finite width a 
and height b , we can sum from x = 0 to x = a the disturbance due to 
a disk located at x = £ whore £ ranges from 0 to a . Notice, how- 
ever, that here one must make a choice of the shape of the loading on 
the blade. That is, we must require that each disk contribute some 
fraction f(§)d§ of the disturbance of the blade. In the problem pro- 
posed in this report, we are at liberty to pick the function £(|) . We 
then go on to find the shape of the blade. However, the inverse problem 
of specifying the blading and then finding the flow is much more diffi- 
cult. In that case, the loading function is in general a function of both 
x and y and depends on the entering flow angles. 

C. Superposition of Actuating Disks for Multi-Stage Compressors 

Given the function f(|) , we can proceed to determine the ve- 
locities u and w throughout the compressor. The function f(£) is 
normalized such that 



For a finite stator of width a , height b , and leading edge at x = 0 , 
the axial velocity for x > 0 becomes 
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For x > a 



u * u + cu - cu e 
o o 




( 22 ) 



0 



The equations for a rotor are similar with only the algebraic sign o£ 
u q changed throughout. 

Equations (20) - (22) give the flow due to one blade row alone. 

For a repeating flow pattern, the velocities are the same as if there 
were an infinite number of blade rows. Thus, for the complete flow pat- 
tern, we must superimpose the flows due to many stators and many ro- 
tors. The following definitions will be used: 

a„ = width of a stator row 
s 

a = width of a rotor row 

T 

c g = gap between stator and rotor 

c = gap between rotor and stator 

r 

f (£)d£ a fraction of disturbance due to stator 

3 fraction of disturbance due to rotor 

£ = a + c + a + c 
s s r r 

The leading edge of the stators are then set at 



x = ml where m = 0, +1, +2, 



9 • • • 



and the leading edgesof the rotors are at 
x = a + c + mi, . 

3 S 

The axial velocity in a stator row due to all blades becomes 
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u = u(y) + cu Q (y) 
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Carrying out the summations, equation (23) becomes 
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(24) 



For the axial velocity in a rotor row, change the sign o£ u q and inter- 
change the s and r subscripts in all terms o£ equation (24). 



-46- 



The corresponding calculation of the whirl velocity in a stator 
results in 

x 
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The whirl velocity in a rotor is determined by changing the sign of u q 
and w and interchanging the s and r subscripts in equation (25). 

An example will be shown later for which the following assump- 
tions are made: 

a s = a r = a * 
c s 3 c r = c • 

ye) » ye> = l/a , 

t = 2(a + c) . 



D. Flow Patterns for Blading of High Aspect Ratio 

If then the exponentials are expanded in polynomials of a/b , 
c/b , or x/b , as the case may be, and terms through second order are 
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retained, the velocities reduce to 

U = U + CU* ( | ^ ) 

— . , x .. ™._x . Trx.ira, 

W = W + ew (2 — - 1) - f-r (2 — - 1 - -r- + -r- ) 
o ' a 7 'a b b 

for a stator, and 

u = 5 - ^u o w { ^ ^ ) 

W = w- rw o (2-- 1) + er (2 - - 1 - ^ + -^ ) 



for a rotor. In both cases 
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The entering and leaving flow angles for both stator and rotor 
blades can now be calculated. The defining equations for these angles 
are, in the case of a stator, 



tan Yj 



u(x = 0) 
w(x = 6) 



and 



tan y £ ~ 



u (x = a) 
w (x = a) 



(29) 



where is the Btator entering angles and stator leaving 

angle. For the rotor, 

tan S = u ( x = 0) 

^ 1 blade speed - w{x=0y 



and 

t ^_ R _ u(* = a) 

2 blade speed - w(x=a) 

In summary, given the blade loading, blade aspect ratio, the 
axial velocity through the compressor, the whirl velocity at one radius 
of a blade row, the desired pressure rise across a stage, and the de- 
sired degree of reaction, one can compute the shape of the blades subject 
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to the restriction that they be of high aspect ratio and very closely 
spaced. 
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m. THEORY FOR FINITE HUB RADIUS 



The analysis for the limiting case of hub ratio equal to one is 
extended to the axisymmetric case of finite hub radius with only slight 
modifications. The y coordinate goes over into the radius r . The 
velocity in the r direction remains v . The major change here is that 
the radial pressure gradients must be accounted for properly. 

The equations of motion in the axisymmetric case are 
Du 1 9 . . n 
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9w , 9w , vw n 

u -5^ + v -3? + — = 0 

As before, we assume a perturbation approach of the form 
u = u(r) 7 eu Q (r) + euj(x, r) 
v = ** v j (x, r ) 

w = w(r) 7 fw Q (r) + rWj(x, r) 
p = p(r) 7 ep Q (r) + RPj(x, r) 



( 32 ) 



Notice that p and p Q are now functions of the radial variable, r . 



The resulting linear equations are 
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and also the radial pressure equilibrium equations 

1 Op ~ 2 
p <5r 



\v 



( 34 ) 
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which must be satisfied at x = + co . 

Eliminating the other dependent variables, the equation for v^ 
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An appropriate solution is of the form 

x 



+ X 



Vj = e 



R(r) . 



Substituting into equation (35) yields 

1 \ 2 1 

R"+.iR' + r (£) - —if - F(r)lR = 0 (36) 

r 

where the prime indicates differentiation with respect to r . If F(r) is 
chosen identically equal to zero, then (36) reduces to Bessel’s equation. 
If we then take only the first term of the solution for v^ , the boundary 
condition that v = 0 at the hub and the tip gives, for the first eigen- 
value 

= 3.18 

2 

for a hub ratio of 0. 6 . If F(r) is small compared to (X ^ /b) , then 
the first eigenvalue does not change substantially. We will now choose 
F(r) identically equal to zero and later compute its actual value to 
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compare with (\j/b)" . 

Proceeding as in the case with infinite hub radius, the perturba- 
tion quantities become 
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where r^ is the hub radius. 

The condition that the total pressure rise across a disk be a 
constant yields 



AP T 

“p 



2c 



‘ o 
P 



+ 2c(uu +ww ) - 2e 
o o 



1 du 
r cTr * 



w / dw 

“ V7I7 

ur 



r / 



riu o (ii)dn.(38) 



Again choosing the special condition of constant total pressure in the 
radial direction implies to first order in c that 
P + P (w + u Z ) = Kj 

and (39) 

p + p (ww + uu ) = . 

1 o r ' o o 2 

Then, as in the simpler case, equation (38) applied across a stator re- 
quires j &2 = 0 . Combining equations (34) and (39) results in the two 



conditions 
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1 u . — . 

-~2 37 ( rw ) 




and 



(40) 




The remainder of the analysis, consisting of the summation of 
the disturbances due to the dicks in order to produce a finite blade and 
the summation over all blade rows to derive the complete flow pattern 
for the compressor, proceeds exactly as in the limiting case, and only 
results will be given here. ‘With the definition 



r , 




the velocity components for the stators are 
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and for the rotors 




and 



(42) 
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The blade angles are determined exactly as before with the use of equa- 
tions (28) and (29). 
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IV. AN APPLICATION OF THE THEORY 



Before proceeding with an example that will allow us to verify 
this theory against an existing compressor, we will indicate the steps 
in order. 

1. Choose the blade loading function f(£) for both the rotor and 
the stator. 



2. Choose the blade dimensions. 

3. Solve equation (36) for the eigenfunctions and determine the 
eigenvalues by application of the boundary conditions on the radial ve- 
locity v . For many applications, F(r) can be taken equal to zero and 
all constants in the series solution for v taken equal to zero except the 
first ones. For such applications, the magnitude of X j will be approxi- 
mately tt . 



4. Continue the summation procedures that result in relations 
similar to equations (41) and (42). 

must satisfy the condition that 

0 
x • 

l 

where r. and r ^ are the inner and outer boundaries of the cylindrical 
annulus. 

6. Determine w from the first of equations (40). There will 
remain an undetermined constant to choose. 

7. Choose the pressure rise across a stage and the degree of re- 
action desired. The pressure rise across a stator 2p Q can then be de- 






5. Pick u and u . Notice u 



, r u dr = 

J o 



ter mined. 
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8. Determine w q from the second of equations (40). The unde- 
termined constant is then evaluated by use of the second of equations (39). 

9. Evaluate the function <$> . 

10. Calculate the axial and whirl velocities for both the stator 
and rotor by using the equations derived in step 4. 

11. Compute the blade angles by use of equations (28) and (29). 

It was decided to use this theory to compare with the blade design 
given in Ref. 1 for solid body blading. The steps will be taken as listed 
above and the results presented here and in Table 1 on page 58. 

1. The loading function was taken to be constant for ease in 
computation. A triangular distribution would be more nearly correct. 

2. The ratio a/b was taken to be 1/3 . 

3. F(r) was taken equal to zero and \ ^ found to be approximately 

3. 18 . The function F(r) was subsequently checked and found negli- 

. 2 

gible compared to (\/b) . 

4. Equations (41) and (42) were then correct as stated in the 

theory. 

5. u and u Q were then chosen to approximate as closely as pos- 
sible the axial flow used in the design of the solid body blading in Ref. 1. 
The relations were 

u/u T = - . 560 | + . 916 

and 

u /u„ = . 185 i - . 151 
o T 

where § is a non-dimensional radius based on the tip radiu3 and u^ is 
the rotor tip speed. 

6. Equation (40) was used to determine w in the form 
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(£) * • 342 1 -. 157 5 2 -!* . 

T | 

The constant Cj was determined by choosing w = . 335 at £ = 0. 8 . 

7. The ideal pressure coefficient for the compressor was 0. 40. 
Then by definition 




0. 40 

~T~ 



u„ 



2 



where u^ is the blade tip speed that has been used throughout this 

analysis in order to non- dimens ionalize the velocity components. Fifty 

per cent reaction in this compressor occurred at £ = 0. 8 . Since 

2 P 0 = P s , we obtained the following relation 

2 

Pq/ p = 0. 05 tirp • 

8. The equation for w q then became 

— = — ! — (. 052 i Z - . 085 % - 2-JJ22 \ . 

U T »/u T ^ IF' 

9. The function $ was evaluated numerically and is listed in 
Table 1. 

10. and 11. The values of u and vv are listed in Table 1 for 
different values of £ . Similarly, the entering and leaving angles y j , 
y 2 , {3j , and are tabulated. The values of these angles used in the 
original design are also tabulated. 
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V. CONCLUDING REMARKS 

For most compressor designs occurring in practice, the aspect 
ratio of the blades is neither very low nor very high. It falls some- 
where in between. The theory of Part I of this thesis is applicable for 
blading of low aspect ratio, where it is assumed that equilibrium flow 
patterns exist at the leading and trailing edges of the blade rows. A 
blade row is not affected by neighboring blade rows except through the 
incident velocity distributions. The theory of Part II for blades of high 
aspect ratio contains the influence of disturbances of potential type that 
can be felt upstream as well as downstream. 

The close agreement between the results of the two theories for 
an aspect ratio of three indicates that there the flow between blades is 
very near its equilibrium state at the leading and trailing edges. 
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